QUANTUM ENERGY EXPECTATION IN PERIODIC 
TIME-DEPENDENT HAMILTONIANS VIA GREEN 

FUNCTIONS 



CESAR R. DE OLIVEIRA AND MARIZA S. SIMSEN 

Abstract. Let Up be the Floquet operator of a time periodic 
hamiltonian H(t). For each positive and discrete observable A 
(which we call a probe energy), we derive a formula for the Laplace 
time average of its expectation value up to time T in terms of its 
eigenvalues and Green functions at the circle of radius e 1 / 7 '. Some 
simple applications are provided which support its usefulness. 
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1. Introduction 

Consider a general periodically driven quantum hamiltonian system 

H{t) = H + V(t) 

with period r acting in a separable Hilbert space 7i and let Up denote 
its Floquet operator, so that if £ is the initial state (at time zero) of 
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the system then U™£, is this state at time mr. Typically, the unper- 
turbed hamiltonian H is assumed to have purely point spectrum so 
that the same is true for e~ tTH °. What happens when H is perturbed 
by V(t)l A natural physical question is if the expectation values of the 
unperturbed energy H remain bounded when V(t) ^ 0. This ques- 
tion is formulated based on many physical models, in particular on the 
Fermi accelerator in which a particle can acquire unbounded energy 
from collisions with a heavy periodically moving wall. Here quantum 
mechanics is considered and, more precisely, if 

sup \(U^,H U^)\ 

is finite or not, where £ € dom H C 7i, the domain of H Q . 

Motivated by models with hamiltonians as above H(t) = H + V(t), 
one is suggested to probe quantum (in) stability through the behavior 
of an "abstract energy operator" which we call a probe operator and 
will be represented by a positive, unbounded, self-adjoint operator A : 
dom A C H — > Ti and with discrete spectrum, 

Aif n = \ n <-Pn, 

< A n < A n+ i, such that if U™£ G dom A for all m E N, then, for 
each m, the expectation value E^{m) = (U™£,, AU™£,) is finite. It is 
convenient to write E^(m) = +oo if U™£, does not belong to dom A. 

We say the system is A- dynamically stable when E^{m) is a bounded 
function of time m, and A— dynamically unstable otherwise (usually we 
say just (un)stable). If the function E^{m) is not bounded one can ask 
about its asymptotic behavior, that is, how does E^{m) behave as m 
goes to infinity? Usually this is a very difficult question and sometimes 
the temporal average of E^{m) is considered, as we will do in this work. 

Quantum systems governed by a time periodic hamiltonian have 
their dynamical stability often characterized in terms of the spectral 
properties of the corresponding Floquet operator. As in the autonomous 
case, the presence of continuous spectrum is a signature of unstable 
quantum systems; this is a rigorous consequence of the famous RAGE 
theorem [T3] , firstly proved for the autonomous case and then for time- 
periodic hamiltonians [121 [25] . At first sight a Floquet operator with 
purely point spectrum would imply stability, but one should be alerted 
by examples with purely point spectrum and dynamically unstable 
[T71 [22l [15] in the autonomous case and, recently, also a time-periodic 
model with energy instability [16] was found. 

Dynamical stability of time-dependent systems was studied, e.g., 
in references [IHl [H [121 [251 I23l [TUl El [H HH [2U [2H [2Ql [2]. In [I] it 
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was proved that the applicability of the KAM method gives a uniform 
bound at the expectation value of the energy for a class of time-periodic 
hamiltonians considered in |18j . 

For hamiltonians H(t) = H + V(t), not necessarily periodic, with H 
a positive self-adjoint operator whose spectrum consists of separated 
bands {(Tj}jLi such that Oj C [Aj,Aj], upper bounds of the type 

(U(t, 0)V>, H U(t, 0)ip) < cte t 2 ^ 

were obtained in [25] if the gaps Xj+x — Aj grow like j a , with a > 0, and 
if V(t) is strongly C n with n > [^r 2 ] + 1. The proof is based on adia- 
batic techniques that require smooth time dependence and therefore do 
not apply to kicked systems. In [231 E] upper bounds complementary 
to those of [25J described above are obtained. 

In [191 [7J H21 E] stability results are obtained through topological 
properties of the orbits £(t) = U(t, 0)f for while in [2Tj |2iH2D[ |2] 

lower bounds for averages of the type 

1 T 

m=l 

are obtained for periodic hamiltonians H(t) = H + V(t) through di- 
mensional properties of the spectral measure /ig associated with Up and 
^ (the exponent 7 depends on the measure /i^). 

In this work we study (in) stability of periodic time-dependent sys- 
tems. As for tight-binding models (see [9] and references therein) we 
consider the Laplace-like average of (Up£,,AUp£), that is, 

m=0 

where A is a probe energy, £ is an element of dom A and Up is the 
Floquet operator. The main technical reason for working with this ex- 
pression for the time average is that it can be written in terms of (see 
Theorem [T]) the eigenvalues of A, i.e., Atpj = Xj(fj, and the matrix el- 
ements (cpj, R z (Up)£) of the resolvent operator R z {Uf) = (Uf — zl)^ 1 
(with z = e~ tE e 1 / T ) with respect to the orthonormal basis {fj} of the 
Hilbert space (here 1 denotes the identity operator). Lemma [T] relates 
the long run of Laplace-like average with the usual Cesaro average. In 
Section [2] we shall prove this abstract results and present some appli- 
cations in Section [3j 

Since our main results are for temporal Laplace averages of expecta- 
tion values of probe energies (see Section [2]), in practice we will think 
of (in)stability in terms of (un)boundedness of such averages. Note 
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that unbounded Laplace averages imply unboundedness of expectation 
values of probe energies themselves. 



2. Average Energy and Green Functions 

Consider a time-dependent hamiltonian H(t) with H(t + r) = H(t) 
for all t G R, acting in the separable Hilbert space H. Suppose the 
existence of the propagators U (t, s), so that the Floquet operator Up = 
U (r, 0) is at our disposal. Let A be a probe energy and Xj, (fj as in the 
Introduction. Also, {(Pj}^i is an orthonormal basis of Ji. 

The main interest is in the study of the expectation values, herein 
defined by 

F A ( ,._S (Un,AU™0, if U^edomA, 
« K ' ' \ +oo, if Uf$ G H \ dom A, 

as function of time m G N. Another quantity of interest is the time 
dependence of the moment of this probe energy, which takes values in 
[0, +oo] and is defined by 

Our first remark is the equivalence of both concepts (under certain 
circumstances). 



Proposition 1. IfUp L £ G dom A for all m, then 

Ef{m) = M/(m), Vm. 

This holds, in particular, if dom A is invariant under the time evolution 
and £ G dom A. 

Proof. Since dom A C dom A^ [13] one has £/™£ G dom A^, for all m, 
and so 

oo oo 

M({m) = x;i<^Vi»^oi 2 =Ei<^i,^>r 

oo 

= Y,\^ Aku F0\ 2 = \\^u^f 

which is the stated result. □ 
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We introduce the temporal Laplace average of (see also the Ap- 
pendix) by the following function of T > 0, which also takes values in 
[0,+oo], 

(1) Lf(T):=|X; c "^H- 

Under certain conditions, the next result shows that the upper ft and 
lower (3~ growth exponents of this average, that is, roughly they are 
the best exponents so that for large T there exist < c\ < c 2 < oo 
with 

ciT p ~ < Lf(T) < c 2 T p+ , 
and the corresponding exponents for the temporal Cesaro average 

m=0 

are closely related; this follows at once by Lemma [TJ which perhaps 
could be improved to get equality also between lower exponents. Note 
that, although not indicated, these exponents depend on the initial 
condition £. 



Lemma 1. If (h(m))™ =0 is a nonnegative sequence, and h(m) < Cm 
for some C > and n > 0, then ft = ft and ft < ft, where 

ft = hmsup MELoM™)) ft = liminf MELoM™)) 

T^oo logT T^oo logT 



ft = hmsup lQ gE T o e -^)) ; (3- = liminf ME^ e-"Mm)) 

Hd T^oo logT ' Hd T^oo logT 

Proof. Note that for < m < T we have e~ 2 < e~^ < 1, and so 

T T oo 

h(m) < ^ e 2 e~~^h{m) < e 2 ^ e~^h{ 

m=0 m=0 m=0 

Hence ft < ft. 



2m 

'ml. 
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On the other hand, for each e > 0, denoting by \x] the smallest 
integer larger or equal to x, one has 

oo nr 1+E l oo 

^^e~~/i(m) = e~~h(m)+ e~~h(m) 

m=0 m=0 m=[T 1 + e ]+l 

\Tl+^ oo 

m=0 m=\T 1 + e ~\+l 

Now, for T large enough ?f < T 1+e < \T 1+e ] . Thus 



OO /»00 

e"— m n < / e-~t n dt 



m =[T i +<i+i " r T1+£ i 

Therefore, for each e > and T large enough 

VV^/i(™) < V/i(m)+C/ e~^t n dt 

|- T l+e] 

< ^ /,( m ) + Ce- 2Te T n . 

m=0 

Since e ~ 2Te T n — > as T — > oo, it follows that 



T^oo lOg 1 

< y logELC 1 Km) 

< lim sup — 

T^oo WgT 

Slin logEltr 1 M"0 logpr 1 -*! 

" TS P logrr^i bgT 



< lim sup 



T ^oo" logfTi+^l logT 

= (l + e)limsup l0gE -=° ClMm) 
(1 + ejlimsup log[Tl+e] 

As e > was arbitrary, f3^ < f3+. □ 

Recall that the Green functions G\{j) associated with the operators 
A, U F at £ G H and zeC, \z\ ^ 1, are defined by the matrices elements 
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of the resolvent operator R z (Uf) = {Uf — zl) 1 along the orthonormal 
basis {ifj}'jL 1 , that is, 

Gl(j) := (ip j: R z (U F )0. 

Note that G\(j) is always well defined since for \z\ 7^ 1 that resolvent 
operator is bounded. Theorem [1] is the main reason for considering 
the temporal averages L^{T). It presents a formula that translates the 
Laplace average of wavepackets at time T into an integral of the Green 
functions over "energies" in the circle of radius e 1 ' 7 in the complex 
plane (centered at the origin). As T grows the integration region ap- 
proaches the unit circle where the spectrum of Uf lives and R z (Uf) 
takes singular values, so that (hopefully) A-(in)stability can be quan- 
titatively detected. 

Theorem 1. Assume that U™^ G dom A for all m > 0. Then 

(2) Lf(T) = r^E A W \ G iO)\ 2 dE, z = e~ iE ^ . 

vre t 1 j=1 Jo 

Before the proof of this theorem, we underline that this formula, that 
is, the expression on the right hand side of (j2J), is a sum of positive terms 
and so it is well defined for all £ G TC if we let it take values in [0, +00]; 
hence, in principle it can happen that this formula is finite even for 
vectors not in the domain of A, where L^(T) = +00. The general 
case, i.e., V£ G TC , can then be gathered in the following inequality 

(3) Lf(T)> r^E A W |GS0')fdE, z = , 

ire t 1 . =1 J 

so that lower bound estimates for this formula always imply lower 
bound estimates for the Laplace average. 

Proof. (Theorem [1]) First note that, by hypothesis, U™f; G dom A? for 
each m G N. Denote by fij the spectral measure of Up associated with 
the pair (ipj, £) and by T the Fourier transform JF : L 2 [0, 27r] — > / 2 (Z). 
By the spectral theorem for unitary operators 

(<Pi,U F $= e- lE 'd N {E'). 
Jo 

For each j let = (a^(m)) m ez be the sequence 

if m < 

e~T J** e - iE ' m dnj(E') if m>0 ' 
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Since E / 1 (Z) fl Z 2 (Z) and T is a unitary operator, it follows that 
= ||^ r_1 aW|| L 2[o,27r] and also 



(jF-y^))^) = — = J2 e lEm a {3 \m) 

V ^TT m=— oo 

= -J—Ye Em e^ / e- iE ' m d N {E') 



m=0 

2tt 



V2^7 l- e *(s-s'+i) ^ 
1 f 2n d^(E') 



1 f 2lT dfijiE') 



1 



27r e lE e t 



-Gl(j), 



with z = e i£+ t. Therefore 



|^- 1 a°')| 2 (E) = -^|GKj)| 2 , 
Zire t 



and so 



/•27T 



IL 2 [0,2tt] 



1 r Z7T 

= - t \Gl(j)\ 2 dE. 

2ne t J 
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From such relation it follows that 



oo „ 



T~ 

m=0 

oo oo _ 

j'=l m=0 

„'_1 J ™_n 1 ^0 



j=l m=0 

00 

VA - lla^ll 2 

^ II 1 IA\ \ \i 



E ll^ 7 a Nl-' 0-2- 



_2_ rxi 



EaJ igiooi 2 ^, 
i=i ■ / ° 



which is exactly the stated result. □ 

Theorem [1] clearly remains true if the eigenvalues Xj of A have fi- 
nite multiplicity. In this case, for each Xj consider the corresponding 
orthonormal eigenvectors cp^ , • • ■ , <fj k , and one obtains 

1 I 00 / k i>2tt n 

Lf(T) = — x^E A J E / \(v^R z (U F )0\ 2 dE 

7Te T j=l \n=l ^° / 

with z as before. 

In case the initial condition is £ = </?i, put t]^- 1 := R z (Up)ipi. Thus, 
(Up — 2)77^ = (pi and so Uprf*' = zrf z ^ + <£>i. Hence 

and by denoting 

G 2 (j):=Gr(j), 

one concludes 



Lemma 2. 

G 2 (j) = 



\ i(^,f/ F ^)), 



if j = 1 
if j>l 
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In Section [3] we discuss some Floquet operators that are known in 
the literature and analyze their Green functions through the equation 



3. Applications 

This section is devoted to some applications of the formula obtained 
in Theorem [TJ In general it is not trivial to get expressions and/or 
bounds for the Green functions of Floquet operators, so one of the main 
goals of the applications that follow are to illustrate how to approach 
the method we have just proposed. 



3.1. Time-Independent Hamiltonians. As a first example and il- 
lustration of the formula proposed in Theorem [U we consider the spe- 
cial case of autonomous hamiltonians. In this case H(t) = H for all t 
and we assume that H is a positive, unbounded, self-adjoint operator 
and with simple discrete spectrum, H {pj = Xj^Pji so that {^j}f =1 is an 
orthonormal basis of TL and < xi < X2 < Xs < ' ' ' with Xj — ► 00 • For 
q > we can consider Hq as our abstract energy operator A, so that its 
eigenvalues are Xj = x] (since A and H have the same eigenfunctions, 
we are justified in using the notation cpj for the eigenfunctions of H ). 
We take U F = e~ iH ° (time t = 1) and for f G H 

Gl(j) = (^R z (H )0 = (itW(tfo)^,0 = 4^-- 

c A ^ — z 

Since dom Hq is invariant under the time evolution e~ ltH °, then for 
z = e~ lE er and £ G dom Hq, we have 

? 1 1 00 f 2n 

L?(T):=Lf°(T) = r^E4 Hti)\ 2 dE 

ire-r 1 Jo 

ire t 1 ~— ^ J \e ^ — z| 

Tie t 1 . =1 j \e v — z| 

Thus we need to calculate the integral L := f„ 2?r , — -r. Let 7 be 
the closed path in C given by j(E) = e lE with < E < 2tt, acj = ere'* 
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and (3j = e re ! %, then 



2 77 



dE 



[e-ix, - z )(e i ^ -z) 
2n dE 



(e"^j - e~ iE eT ! )(e i ^ - e iE eT) 
dE 

2 i '■ '. I '■ '. 

o er(e~?e"'* — e~ lE )(e~T e% — e lE ) 



2tt 



2t 



dE 



Jo e- iB e~^e-^ (e iB - aj)(e iE - (3j) 



2- 



ie iE dE 



ere %x 3 i Jo 
i 



{e iE _ a . ){e iE _ Pj) 

dw 



e^e-vo J 1 (w - aj)(w - fy) ' 
As \aij\ > 1 and \(3j\ < 1, /3j is the unique pole in the interior of 



Thus, by using residues, 



-2m- 



1 



2ty 



eTe-^ {Pj-atj) e f _ i 



and Ij is independent of Xj- 
Therefore by (jlj) it follows that 



L%T) 



2_ rri 

ire t i 



2 1 



e tT ( e i _ 1 



1 oo 

--rE>JI<yi,Ol s 



1 — e r 



T 



Since ^1 — e rj = | ; + 0(J^) ; f or large T it is found that 



L\{T) 



with (for £ G dom Hq) 



Then we conclude that the function 



lim L\{T) = ^HU). 

— >OQ 

he function 
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is bounded for £ e dom Hq, which is (of course) an expected result (see 
Proposition [Q. 



3.2. Lower Bounded Green Functions. As a first theoretical ap- 
plication we get dynamical instability from some lower bounds of the 
Green functions. See [9] for a similar result in the one- dimensional 
discrete Schrodinger operators context; there, a relation to transfer 
matrices allows interesting applications to nontrivial models, what is 
not available in the unitary setting yet (and it constitutes of an impor- 
tant open problem). As before, Xj denote the increasing sequence of 
positive eigenvalues of the abstract energy operator A, the ones we use 
to probe (in) stability. 

Let [•] denotes the integer part of a real number and | ■ | indicates 
Lebesgue measure. 



Theorem 2. Suppose that there exist K > and a > such that for 
each 2N > large enough there exists a nonempty Borel set J(N) C S 1 
such that 



\Gi(j)\ N<j<2N, 



N- 



holds for all z = e~ iE+ ^ with E £ J T (N) = {E" 6 S 1 : 3 E' G 
J{N); \E" - E'\ < J} (the ^-neighborhood de J(N) ). Let 5 > 0; then 
for T large enough such that N(T) = [T 5 ], one has 



Lf(T) > cte Ap^ 1 - 2 ^- 2 . 
Moreover, if Xj > cte p , 7 > 0, then 



Lf{T) > cte T 5 ^~ 2a+1) - 2 . 
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Proof. By the formula in Theorem (TJ or its more general version 

OO „ 2 7r 



1 1 J— f 

Lf(T) > r^E A ^ / \ Gi M dE 

ire t i Jo 



j=JV(T) 

27V(T) 

> £ / |Gla)| 2 dE 

J j= 7V(T) J ^( JV ) 

2N(T) 

cte , v-^ ^ i T 

> -^co £ ^^| Jt (jv)| 

j=N(T) K 1 

cte | , * | K 2 

| Jt{^ )\^N(T) 



= ^-|J r (iV)|A [T5] a _ x 

> cte A^T^ 1 - 2 ")- 2 ; 

we have used that | Jx(iV)| > ^. If Aj > cte p then 

Lf(T) > cte r^T^ 1 - 2 ")- 2 = cte T s ^ 2a+l ^ 2 . 

The proof is complete. □ 

The above theorem becomes appealing when the exponent of T is 
greater than zero and instability is obtained, for instance when 5(7 — 
2a + 1) > 2 in case Aj > cte p . However, up to now we have not yet 
been able to find explicit estimates in models of interest; in any event, 
we think it will be useful the future applications and so we point out 
some speculations. First, note that it applies even if the set J(N) is a 
single point! Nevertheless, we expect that Theorem [2] will be applied to 
models whose Floquet operators have some kind of "fractal spectrum" 
(usually singular continuous or uniformly Holder continuous spectral 
measures) and, somehow, a should be related to dimensional properties 
of those spectra; indeed, this was our first motivation for the derivation 
of this result and, in our opinion, such applications are among the most 
interesting open problems left here. 
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3.3. Rank-One Kicked Perturbations. Now consider 
H(t) = H + kP^ 5 (t ~ ™ 27r )> 

n 

with H as in Subsection I3.1[ with eigenvectors {^fj}f =1 and Xj the 
corresponding eigenvalues; P^(-) = (<j>, -)4> where k G M and is a 
normalized cyclic vector for H , in the sense that ||0|| = 1 and the 
closed subspace spanned by {H™(j) : m 6 N} equals H . Let 

3 

In this case (see [El E]) 

U F = U (1 + «P^) , 

with C/ = e - j27r ^° and a = (e" i27rft -l). Note that <\> e dom ifj, Vg > 0, 
and so for £ e dom ifj, 

also belongs to dom Hq] a simple iteration process shows that U™^ e 
dom Hq for all m > and we are justified in using the formula in 
Theorem [1] to estimate Laplace averages. 

We are interested in rf z > = R z {Uf)<£i- As \z\ ^ 1 it follows that rf x ^ 
belongs to the Hilbert space and so one can write 

oo 

Note that aj = G z (j) and we have 

(5) tW* - ^ W = Vi- 

By the relation 

oo 

= ^ajC/o^j + a[/ o (0,r/ (z) )(/) 

oo oo 

= E + «<0, E v ,2;:v v, 

j=i j=i 
= E(% + afav^e-^ips, 
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and (J5]) it follows that 

oo oo 
3=1 3=1 

that is, 

oo 

He~^ -z) + a^d z) )^e~ l2 ^\ Vj = Vi, 

3=1 

and we get the equations 

ai ( e -i2*Xi _ ^ + ^ ^W^g-iftrxi = 1; 

aj{e- i2 " x * -z) +a(<f),r] {z) )b J e- i2n ^ = for j > 1. 

Thus 

1 - aifarj^he-* 2 ** 1 



(6) ai 



e -i27r X i _ ^ 



a(0,7yW)6 J -e- i2 ^ . 

For the trivial case a = or, equivalently, kgZ, one has 

1 



e -i2vrxi _ ^ ' 

% = 0, J > 1, 
and rf z > = e -atx 1 _ z - In this case the analysis of L^(T) is reduced to 

2 " 2 .„ f 2n dE 2tt 



o ■ ■ Jo |e-^ Xl - ^ 2 - 1 

as calculate in Subsection 13.11 Thus L*(T) ~ Hf/oVill f° r large T, as 
expected. 

Returning to the general case a^O, note that 

oo 



- / I -q(^,7y (2) )6ie- i2,rXl \ y,- (-q) (0, rj^bje-' 2 ** 

j=2 



bi a\b-\ 2 e~ i2nX: > 

3=1 
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So 

a\bj\ 2 e- i2lTX ^-\ 

(r \ ' - : ) 7 



• ^ = ( P -i2*Xi - z) 1 + ^ 



By denoting 

T ( z ) = 1 + v 1 jl 



by ([6]) and ([7]) we finally obtain the relations 



1 a 

CL\ = 



h x \ 2 e~ l2 ^T(z)- x 



a6 j 6ie- i2 ^r(^)- 1 . 

3.3.1. A Harmonic Oscillator. Now we present an application of the 
above relations to a kicked harmonic oscillator with natural frequency 
equals to 1; we will write L q ^ = . 



Proposition 2. Let H Q be a harmonic oscillator hamiltonian with ap- 
propriate parameters so that its eigenvalues are integers j, j > 1, and 
Uf = Uq(1 + aPt/)) as above. Then for any kg! and cyclic vector 
for H , there exists C > so that, for T large enough, 

L q vl (T)<C, 

where tp\ is the harmonic oscillator ground state. Hence we have HQ- 
dynamical stability. 

Proof. We use the above notation; note that <p\ G dom Hq, Wq > and 
Theorem [I] can be applied. In this case we have 



a|6j| 2 _ a 9 1 — z + a 



and so 



1 a 1 61 

a 1 



1-z (1 - z)(e- i2nK - z) 
abjbi 

a i = ~ (l - z ){e- i2 ™ - z)' 3 > 
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2tt 



Now we evaluate Ij := f Q n \aj\ 2 dE. For j > 1 and 7(-E) = e ijB , < 
E < 2vr, 



p2w 
JO 



\aj\ 2 dE = 



2n 



abjbi 



(1 - ^)(e- i27rK - z) 

(■2-K 

a\'%\'% 



dE 



dE 



(1 - e - iE eT)(e- i2nK - e- iE er^ 
M 2 |kj| 2 |^i| 2 f w dw 



Jj(w- ft)(w - 



ft)(™-ft)(w-ft)' 

where ft = e^, ft = e~^, ft = e ^e* 27m and ft = e ~^e i27rK ; only (3 2 
and /3 4 are poles in the interior of 7. By residue, for j > 1, 



4 = 



27r|a| 2 |6 j | 2 |6 1 

g g 



x 



+ 



ft 



(P2 - ft) (ft - ft) (ft - A) (A - A) (ft - ft) (ft - ft) 
2Tra\b j \ 2 \b 



2|7T|2 

1 



2 7 ra|& i | 2 |& 1 |V 2 ™ 



(e.T — l)(e * 2?rK - er) (ef - l)(e l27rK — er) 



27ra|6 7 f|6 1 12 



rD i2'KK 



and for j = 1 



/1 = 



2tt 



a|6i 



(l-z)(e-< 



%2-KK 



n2n 
JO 



dE 



o (l-z)(l-z) 



— a; 1 61 



2tt 



(l-z)(l-z)(e i2 ™-z) 



— a\b 



2n dE 
11 7 (l-z)(l-z)(e- i2 ™-z) 

dE 



p2n 

+ l«l>i| 4 / 

Jo 



(1- z)(l- z)(e~ i27rK - z)(e i27TK - z) ' 
evaluating the integrals we obtain 
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and after inserting this in the expression of the average energy we get 

a\bi\ 2 



L U T ) = — zr- i-N - 



2 



2|&!| 2 

2a\b 1 \ 2 / 1 e i2 ™ 



— f v 7 ; ^()( 

Therefore, for large T there is a constant C(k, b\) > so that 

This completes the proof. □ 

For harmonic oscillators with eigenvalues uuj, u ^ 1, the evaluations 
of the resulting integrals are more intricate and were not carried out. 



3.4. Kicked Perturbations by a V in L 2 (5' 1 ). 

3.4.1. Kicked Linear Rotor. Consider 

H(t) = ujp + V(x) S(t - n2n), 

where p = -i£, uj G R and V G L 2 ^ 1 ). The Hilbert space is L 2 (5 1 ); 
this model was considered in [4J [101 [11] and references therein. The 
Floquet operator is 

U F = U V = e -^v e -iV(x)_ 

Denote by (fj(x) = e^ x /\rM, < x < 2n and j £ Z, be the eigen- 
vectors of p 2 whose eigenvalues are the square of integers j 2 ; all eigen- 
values have multiplicity 2 (the corresponding eigenvectors are ipj and 
ip-j), except for the null eigenvalue which is simple. 
Consider the case u = 1; then 



({u F - *)-v„) (x) = 




and so 
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-2-k 



Denote Ij : = f Q \G^°(j)\ 2 dE. It follows that 



2tt 



2-7T —ijx 

dx 



-iV(x) _ z 



(2tt) 2 i 7 Vio (e-^)-z)(e^)-z); 

d£ 

( e -iV(x)_ z )( e jV(y)_ 2 ) 



For x^eS 1 fixed denote := J 27 " (e -iv W _ff (e ^(«)-^) - If = 



< < 2tc, one has 

" 27r dE 



Ixy 



I 



( e -iV{x) _ e -iE £ j!^ e iV(y) _ g«£gT) 

f 27T dE 

JO e -iE e -iV(x)^ e iE _ e iV(x) e ^ ^g^r ^ giV (y) _ e iE^ 

dw 



and by residues 

2tt 2tt 



Hence 



/. — / / e -ijx e ijy dxdv 

3 (2tt) 2 J J ( e f - e -iV(x) e iV(»)) 

(8) = J_ e -0-x ( ^ >) ^ 

27r7 V7 (g| -e-iVWgiVCy)); 



1 r 2 ^ g-ijas , /-27T e ijy dy 



dx. 



2-n Jo e~ iV ^\Jo ( e Te iV W - e iV ^)' 

The analytical evaluation of these integrals is not a simple task. As 
an illustration, consider the particular potential V(x) = x; since by 
Cauchy's integral formula 

f 2n e l3y dy 1 f wi-^dw n . ^ , 

/ T^T" ^ = ~~/7 ^^ = ' lf 3 > 1, 

Jo (ere™ — e l?/ ) 11 Jj [w — ere™) 

and by residue theorem 



Z" 2 *" e ij '"dy 1 f dw 2tt 

it is found that 



er e lx — e lv ) l J^w l i(w — eTe lx ) (ere 81 ) 1 •? 

= if j > 1 



, if J < 0, 
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and 

1 f2* e -ij* 2 7T , 1 f 2lT , 2tt 



Ia — — / — ; = dx = — o- — r / dx — — ~- — -, if 7 < 0. 

Therefore, by ([3]) it follows that for any q > 

and we conclude that (see the Appendix) 

^ 2 9 ( m ) ^ cte m2q 

and also that the sequence m \— > {U^ipQ^p^Up^ipo) is unbounded. This 
behavior is expected since the spectrum of Up is absolutely continuous 
in this case [2], but here we got the result explicitly without passing 
through spectral arguments, although in a rather involved way; indeed, 
a much simpler derivation is possible by direct calculating U™<po and 
the corresponding expectation values. 

For V(x) = kx with integer k > 2, similar results are obtained, that 

is 

if j = Ik, I > 1 



2- 



e 2 / T (l-l) 



if j = lk, I < 



and so 



Z=l 



T 

Therefore, we have the following lower bound for the Laplace average 

L^(m)>C(k,q)m^ 

(see Appendix). The same is valid if V(x) = kx with k denoting any 
negative integer number. 

3.4.2. Power Kicked Systems. Due to the difficulty in evaluating the 
integrals in (JSj) , in order to estimate (T) in some situations we take 
an alternative way. 

Consider the Kicked models in L 2 (S' 1 ) with Floquet operator 

(9) U F = U v = e -v™f(p) e -iV(x)^ 

corresponding to the hamiltonian 

H(t)=u J f( P ) + V(x)Y J Kt-^n), 
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with p, V, tpj as before and f(p) = p N for some N G N. Let T : 
L 2 ^ 1 ) -> / 2 (Z) be the Fourier transform. Then J=\J V J=- X : Z 2 (Z) -> 
Z 2 (Z) and 

TVyT~ x = J^e' 12 ^^^ e~ iV ^ T^ 1 = j^e~ i2nuj f^J-'~ 1 J-'e~ iV( - x " ! J-'~ 1 

where j?e- l2 ™f(j>)jr-i is represented by a diagonal matrix D whose 
elements are 

D(m,n)=e- i ** uf W6 mn , 
and Te~ tV ^T~ x is represented by a matrix W whose elements are 

W(m, n) = (J-'p) (m — n) = p(m — n) , 
where p(x) = -}=e~ lV ( x > . Denote B = DW; so 

B(m } n) = e - i2 ™ f{n) p(m - n) 

and 

(10) U v = T~ X BT. 

Put rf z) = R z (U v )(fo; then 

U vV (*) _ zr f») = ^ 

and using ffTUl) we obtain 

BTii z) - zTii z) = Fif . 

Thus, for each iigZ, 

(BF V W)(n) - (zF V W)(n) = (F<p )(n), 

so that 

(11) e -^fMJ2p(n-j)Gt»(j) - zGT(n) = 5 n0 . 

Tridiagonal Case In order to deal with the above equations, we try 
to simplify them by supposing that V is such that p(m — n) = if 
\m — n | > 1. Then, for each n G Z fixed (TTTj) becomes 

(12) e -«™/(») £ p{n-j)G?{j)-zG?{n) = 6 vD 

\n-j\<l 
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and T x UyT = B is tridiagonal and has the structure 



B 



\ 



-1)P(0) 



(-i)pH 
p(o) 



g(l)p(0) g(l)p(-l) 
g(2)p(l) g(2)p(0) 



\ 



J 



where g(n) = e - i27TUjf( - n \ 

Now, a tridiagonal unitary operator U on / 2 (Z) is either unitarily 
equivalent to a (bilateral) shift operator, or it is an infinite direct sum 
of 2 x 2 and lxl unitary matrices, as shown in Lemma 3.1 of [6]. 
For proving this result it was only used that U is unitary and Uek = 
a^efc-i + j3k£k +7fcefc+i, where {e*;} is the canonical basis of / 2 (Z), that 
is, 



/ 



U 



\ 



ft 
7fe 



7fc+i 



It then follows that for all k G Z 

Kl 2 + |ft| 2 + |7fc| 
7fc_i/3 fc -i + ftSfc 



0. 



1. 



Applying these relations to B — T~ l UvT we obtain 

• If p(-l) ^ then p(l) = p(0) = and |p(-l) 

• If p(l) + then p(-l) = p(0) = and |p(l)| = 

• If p(0) ^ then p(l) = p(-l) = and |p(0)| = 

The next step is to investigate these cases. If p(0) ^ it reduces to 
the autonomous case H(t) = H previously considered. 

The cases p(— 1) ^ and p(l) ^ are similar, so we only discuss 
that p(l) 0. For n£Z fixed, equation (fT2|) takes the form 



1. 
1. 



(13) 



-i2irujf(n) 



p(l)Gt°(n-l)-zG 



Pa , 



n 



6, 



so we can write G^°(n) in terms of G^°(0) and G%°( 
More precisely 



-1) for all n G Z. 



-i2™(/(n)+-+/(l)) 
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~n— 1 

2 v y e -t27rw(/(-n+l)H h/(— 1); p(l ) n_1 — 

moreover, for w = in (jI3j| we obtain p(l)G£°(-l) - zGf (0) = 1, so 
for z = e~ iE e 1/T and T > 1 

1 < \Gt»(-l)\ + \z\\GT(0)\ 
= |Gf(-l)l + e 1/T |Gr(0)| 
< e (|GT(-l)l + |GT(0)|), 
and there exists d > so that 

|Gf(-l)| 2 + |G£°(0)| > d> 0. 
Therefore, by d3J), for T > 1 one has 



1 1 . , / /" 27r /" 27r \ 

> r^J> 2 M/ \GT(n)\ 2 dE+ \G^(-n)\ 2 dE) 

00 / i /.2tt 

VeT Jo 



1 1 



_2_ rxi 

ire t i 



71=1 

2tt 



-I -I 00 /-27T 

^ — (iGr(o)i 2 + |Gr(-i)i 2 )^ 

7re t J n=1 y 

n=0 

so that, by the discussion at the end of the Appendix, 



l£(m)>C{m + l 



i 2( / 



and (Uy(po,p 2q Uy(po) is unbounded. Hence we have instability. 

Pentadiagonal Case Suppose now that V is such that p(m — n) = 
if \m — n\ > 2. Then for each n G Z fixed, equation ffTTj) becomes 



(14) c -«™/(») ^ p(n-j)Cr(i)-zGr(n) = ^o, 

|n-i[<2 

and T~ x UyT is pentadiagonal and has a structure similar to the corre- 
sponding operator in the previous case, just adding the elements whose 
distance to the diagonal is 2. The elements in the new upper diagonal 
are e~ l2nul ^ n ^ p(—2), and the new lower diagonal are e - t2nuJ f( n ) p(2). 
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For not repeating the tridiagonal case we suppose that either p(2) or 
p{— 2) is different from zero. If U is a pentadiagonal unitary operator 
in Z 2 (Z), that is, Ue k = ( k e k _ 2 + a^k-i + Pke k + 7fc e fc+i + d k e k +2, one 
gets the matrix representation 



U 



\ 



Pk-2 
lk-2 
@k-2 



CKfc-1 
Pk-1 

7fc-i 
Ok-i 



a k 
Pk 
Ik 



lk+1 



Ck+2 
fik+2 



\ 



J 



From this we obtain the following relations, for each k G Z, 
|Cfc| 2 + |afc| 2 + l/5fc| 2 + l7fc| 2 + W = l 



Cfc«fc-i + + Pkik-i + ikOk-i = o 



P k -i9 k -i + a fc 7 fc + Ck+iPk+i = 



a k -i9 k ^i + Cfc7fc = 



Suppose that p(2) 7^ 0. The case p{— 2) 7^ is similar. Then by the 
above relations we obtain p(— 2) = p(— 1) = p(0) = p(l) = 0, and so 
( TH|) becomes 

e - i27r£j/(n) p(2)Gf (n - 2) - zGf (n) = <W 

For n = one gets p(2)G£°(-2) - zG£°(0) = 1 and analogously to the 
previous case 

|Gf(-2)| 2 + |Gf(0)| 2 >d>0, 
with z = e~ tE e l l T and T > 1. Since for n > 1 



Gf(2n) 



-i27ru,(/(2n)+/(2n-2)+-"+/(2))^2) n 



Gr(o), 



Gf(-2n) 



z n-l G ^o(_ 2 ) 



p(2) n - 1 e- i27rw (- f (- 2 ( n - 1 ))+-+- f (- 2 )) ' 
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we obtain 

l£(T) > —T-ypn^ll \GT(2n)\*dE+ \ |Gj»(-2n)| a dE 



1 1 


2 

7re t 


T 


1 


1 


2 


T 


2(n- 


-i) 


+e i 




1 


1 



2vr 

|G?°(-2)| 2 dE 



OO ^TT 
2n 



^ T " n=l 



> — yE(2n)VT / (|Gr(0)| 2 + |Gr(-2)| 2 )^ 



> d-£(2(n + l)) 2 V 

hence 



2n 
T 



m 

n=0 



iZ 2 o 9 (T)>C(2(m+l)) 



2g 



and (Uyipo,p 2q Uy(po) is unbounded. 



iV-diagonal Case If V satisfies pirn — n) = for |m — n\ > N, we 
suppose that either p(N) or p(—N) is different from zero. In case 
p(N) 7^ 0, by unitarity and the structure of T~ x XJyT we obtain that 
p(N - 1) = • ■ ■ p(0) = p(-l) = • ■ ■ = p(-iV) = 0, thus dn} becomes, 
for each nGZ, 

e -«f~/(«) i 3(jV)G^ ) (n - JV) - zG£°(n) = 5 n0 , 

and so 

|Gf(-iV)| 2 + |Gf(0)| 2 >d>0, 
with z = e~ lE e l ' T , T > 1. Moreover, for n > 1 

GTinN) = ^^Gf (0) 

and 

Z n-1QW(_N) 

Gf n (-n.\) ■ 



p( i\T) ?1 - 1 e- i27rw W(- Ar ( n - 1 )) + - + ^(- 7V )) ' 
Similarly to the previous cases we conclude that 



o 00 



T 

n=0 



Therefore we can stated the following result: 
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Theorem 3. For Kicked systems in L 2 (5' 1 ) with 

JJ — e -i2™f(p) e -iV(x) 

as in (TjJ), we obtain that TUyT~ x : / 2 (Z) — > / 2 (Z) is represented by the 
matrix B with elements B(m,n) = e~ l2nul f( n ^ p(m — n), where p(x) = 
{2ix)~2e~ lV ( x \ If V satisfies p(m — n) = for \m — n\ > N G N* and 
either p(N) or p(—N) is different from zero, then V(x) = ±Nx + 9, 
for some and TUyT~ x is unitarily equivalent to T N (the Nth 

power ofT) where T is the bilateral shift. Furthermore, 

9 oo 

<(T)>4E(^( n+1 )) a9e ""- 

n=0 

Proof. It is enough to prove that TUyT~ x is unitarily equivalent to 
T N . Suppose that p(N) ^ (the case for p(—N) ^ is similar); then 
by the above discussion we obtain 

, v f ii m^n + N 

^^-\ e -^m m if m = n + N > 

that is, Be n = e -* 27ra; /( n ) p(N)e n+ N where {e n } is the canonical basis 
of 1 2 {Z). Since \p(N)\ = 1, write p(iV) = e~ ie . Let W be the unitary 
operator defined by 

We n = e^ n e n , n G Z, 
where i? n are elements in [0, 2ir). If $ n satisfies for all n G Z 
(15) tf„ +JV - #n = 2vrw/(n) + 0, 

it follows that W~ l BW = T . (II 5p is satisfied taking, for example, 
i?o = i?i = • • • = i?jv_i =0 and the another obeying (IT51) . □ 

Although Theorem [3] gives a nice illustration of the potential appli- 
cations of our expression for the Laplace average, since it is one of the 
few instances that such average can be explicitly estimated from be- 
low, again it can be derived by more direct methods and one can also 
conclude [3] that the spectrum of the corresponding Floquet operators 
are absolutely continuous. 



4. Conclusions 

Although most of our applications of Theorem [1] give expected results 
(sometimes known results that can be derived in simpler ways), we 
believe that that formula is interesting and has a potential to be applied 
to more sophisticated models as the Fermi accelerator. The difficulty is 
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to get expressions or estimates for the Green functions, since calculating 
the resolvent of an operator is not always an easy task; sometimes we 
have the expressions for resolvent operators (e.g., for kicked systems) 
but the resulting integrals can be too involved. We have not tried any 
numerical approach to formula (j2J), which might be useful for some 
specific models. 

In the case of one-dimensional discrete Schrodinger operators, where 
the hamiltonian is Hy : / 2 (Z) — > Z 2 (Z) defined by 

(H v Z){n) = £(n + 1) + £(n - 1) + V(n)£(n), 

V a bounded sequence, a similar formula can be handled in some cases 
by relating the resolvent R E+ ±{Hv) to transfer matrices. Then ade- 
quate upper bounds of such transfer matrices, on some set of energies 
E, result in lower estimates for the corresponding Green functions and 
then transport properties are obtained for interesting models (see [9] 
and references therein). 

In [6] a class of Floquet operators displaying a pentadiagonal struc- 
ture was introduced; for these models there is a transfer matrix formal- 
ism. However, such transfer matrices are too complicated and analyti- 
cal estimates seem far from trivial. 

Anyway, the technique here is quite general, it asks no particular 
regularity of the time-dependence and can be virtually applied to any 
time-periodic system as soon as the time evolution is well posed. As 
already said, the chief difficulty is related to suitable bounds of matrix 
elements of the resolvents of unitary (Floquet) operators, a task harder 
than we initially envisaged. Herein we put forward for consideration the 
challenge of getting additional applications for the formula fl2]) deduced 
for the Laplace averages, including an application of Theorem [2] to 
physical models. It is also worth mentioning the question left open in 
Lemma [TJ that is, is it true that j3~ = 



Appendix: Laplace Transform of Sequences 

Let a = (a n ) ng N be a sequence of positive real numbers. The Laplace 
transform of a, denoted by f a , is the function defined by 

oo 

(16) /.OO =!>-*» a(n), 

n=0 

for s in a subset of E. It will also be denoted by f a {s) = C(a). 
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We say the Laplace transform of a = (o n ) exists if the series in 
converges for some s. For example, if a(n) = e n , then the sum in (fl6l) 
diverges for all s G R. 

Examples. 

(1) For the constant sequence a{n) = 1 it follows that 

oo 1 
sn _ 



n=0 



(1-e- 



for s > 0. By using Taylor expansion, for small s one finds that 
i ~ i 

(l-e- s ) ~ s> ' 

(2) Since ^]^L £ n = tzt, for 2; G C, \z\ < 1, it follows that 
X)(n + fc)(n + fc-l).-(n + l)z"= (1 _ z)k+1 , 



n=0 



for A; = 1, 2, 3, • • • , and z as above. Thus, the Laplace transform 
of a k (n) = (n + k)(n + k — 1) • • • (n + 1) is 

°° u 
f k a {s) = Y j e- sn a\n) = I a > 0. 

n=0 ^ ' 



For small s, /„(s) 



A sequence of complex numbers a = (a n ) is said to be exponential 
of order do (real) if there exists M > so that |a(n)| < Me a ° n , Vn. 
That is, a(n) does not increase faster than e a ° n as n — > 00. If a = (a n ) 
is exponential of order cr > 0, then 

00 

n=0 

is convergent for any s > <Jq. 

Let V denote the set of positive sequences of exponential order a . 
The Laplace transform C satisfies 

C(ca) = cC(a), C(a + b) = C(a) + C(b), 

where c is a positive number and a and b are sequences in V. Moreover, 
if a G V and £(a) = 0, then Y^=o e~ sn a(n) = and so a{n) = for all 
n, that is, a = 0. Thus £ is injective on V. 
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The Laplace average ([I]) is related to the Laplace transform of (n) 

by 

j*m = !E e -»^(») = 

n=0 ^ ' 

If a(n) = 1 for all n, then 



T Ja \T) T(l-e- 2 / T ) ~ T2/T 
for T large enough. If a(n) = (n + fc)(n + k — 1) • • • (n + 1) ~ n fc , then 

T /a v T; T(l - e -2/T)fc+i ~ T(2/T) fc + 1 ~ V 2 / ' 

for large T. Hence, if E^(T) grows like T fc then the same law holds for 
its average Laplace transform. We have a restricted converse, that is, 
if L^{n) grows with a positive power of n then, by Lemma [H its Cesaro 
average is unbounded (with a rather similar behavior at large times) 
and so is E^{n). These properties are repeated used in the text. 

One should be aware that there are special situations of unbounded 
positive sequences a(n) with bounded average Laplace transforms (so 
that /3+ = = 0); an explicit example is a(n 2 ) = n and a(n) = 
for n £ {k 2 : k G N}. The same phenomenon is well known for Cesaro 
averages and, by Lemma [U such phenomena are connected. 
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